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ABSTRACT 



The effective action induced by fermions in the chiral Schwinger model with charges (3,4,5) is 
investigated. Pauli-Villars regularization is combined with momentum cut-off for the evaluation 
of the fermion determinants on continuum gauge fields interpolated between lattice points. The 
convergence and gauge variance are studied numerically on gauge configurations taken from 
quenched updating. 
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1 Introduction 



The non-perturbative definition of chiral gauge theories is a long-standing problem of quantum 
field theory on the lattice (for reviews see, for instance, [1] -0]). Recently there have been 
several interesting new proposals trying to circumvent the obstacles represented by the Nielsen- 
Ninomiya theorem || in different ways Hj— pf (for further references see these papers). 

In the present letter the idea put forward by 't Hooft || is applied to the simple test case 
of the two-dimensional chiral Schwinger model. It is based on the interpolation of the lattice 
gauge field and the definition of the chiral fermion determinant on the obtained continuum 
gauge field, by exploiting the knowledge accumulated in continuum approaches (for a review 
see |JL0[|). Similar ways of defining chiral gauge theories on the lattice were discussed for some 



time [ TT| , p~2[] , and have been recently further developed in refs. [[13|, [14 



The aim of the present paper is to study numerically the definition of the effective ac- 
tion induced by chiral fermions on the interpolated lattice gauge field. As a first step, the 
two-dimensional massless chiral Schwinger model is considered here. This and related two- 
dimensional models are often used as a testing ground for chiral fermions (for recent examples 
see ||15||). Since this model is well known and exacly soluble (in the extensive literature see, 
for instance, the papers in ref. |16[]), the questions are mainly oriented towards the qualitative 
behaviour of the calculation of the effective action along the line of refs. [Q, [T^, H|. The meth- 



ods used will be such that they can be extended to four dimensions in a straightforward way. 
In the next section the interpolation of the U(l) gauge field is discussed. This is followed by 
a short discussion of some useful numerical algorithms. In section [| the convergence of the 
chiral fermion determinant is considered by removing the momentum cut-off. The variation 
with respect to gauge transformations is numerically investigated and discussed. 



2 Gauge field interpolation 

The lattice gauge field is defined by the parallel transporters on the discrete links of the lattice, 
which is chosen in the present paper, for simplicity, to be hypercubic with periodic boundary 
conditions. To extend the gauge field connection into the meshes of the lattice is highly ar- 
bitrary. In order to reduce this arbitrariness some guiding principles must be respected, such 
as smoothness and some minimality principle which chooses among the different possibilities. 
Since in chiral gauge theories the anomaly plays an important role, one can also connect the 



gauge field interpolation to the geometrical definition of the topological charge [|1y], |T^]. Both 
this interpolation and the piece-wise linear one minimizing the Euclidean action, which has 
been proposed in M, have the important property that in momentum space the support of 
the Fourier transform is concentrated near the momenta allowed for the gauge field on the 
lattice. This has to be required as a condition for any reasonable interpolation: the momentum 
cut-off imposed on the gauge field by the lattice has to be approximately maintained by the 
interpolated gauge field too. 

Let us denote the U(l) gauge link variables on the lattice by 

U Xfl = exp(igA Xfl ) . (1) 
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Here x denotes lattice points: x = (xi,X2) with integer x^ satisfying < x^ < L^ — l, (// = 1, 2). 
The lattice extensions are denoted by L M , g is the bare gauge coupling and the number of lattice 
points will be denoted by Q = L X L 2 . Note that throughout this paper the lattice spacing of 
the lattice for the gauge field is set to be a = 1. In other words, every dimensional quantity, 
as for instance the gauge field A x/X , is measured in lattice units of the gauge field lattice. The 
Fourier transformation to momentum space is defined, as usual, by 

\ = 5;e^% . (2) 

X 

The inverse relation is ^ 

A-xu, = tz e lk x A kVh , (3) 

11 k 

where the sum is running, of course, on the points of the Brillouin zone 

K = T-^ > ^ = -mt(L M /2), -int(V2) + 1, ... , +int(L M /2) . (4) 



The above discussed condition on the interpolation means that, if the Fourier transformation 
of the interpolated gauge field is performed on the continuous torus, the Fourier coefficients are 
approximately the same as in (0). This suggests the introduction of eq. (|3|) as the definition 
of the interpolation by simply extending its validity to continuous x. This means that the 
interpolated gauge field on the continuum function of the continuous x r \ 

Afj,(x c ) = J2D^(x c - x)A Xfl (5) 

x 

with the interpolation kernel 

D„(x c -x) = ^-J2 eiHXc ' x) ^ kfl ■ (6) 

k 

This is a very smooth interpolation indeed, since the result is infinitely many times differentiable 
(an entire function for complex x c ) satisfying, for integer x: 

A^x + jl/2) = A x ^ . (7) 

Here, as usual, fi denotes the unit vector in direction jj,. 

Before going further, let us make a short technical remark. The Fourier coefficients in eq. (0) 
have the periodicity properties 

Ak+2-Kfi,^ = ~~ Ak,fi , A k+ 2iTu,^ = Ak^ , (8) 

where v = 3 — fi. This allows us, for instance, to choose the symmetric interval for momenta in 
(|J). Of course, for the summation in the definition of in (|6]) one has to specify the interval. 
In fact, for even L M it is advantageous to stick to an exactly symmetric definition by dividing 
the Fourier coefficient at = into equal halfs at = and = —L^/2. (For an 

odd L M the points in (||) are automatically symmetric.) In the numerical study discussed in 
section |] this symmetric definition was always taken. 
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It must be emphasized that the interpolation defined by eqs. ([5]) and (|6]) is only one example 
among many others. It satisfies the important condition that the momentum cut-off for the 
gauge field be transferred from the lattice to the continuum. This is the only condition which 
will be exploited in what follows. This implies that the conclusions from the numerical study 
in section 01 will be qualitatively valid also for the interpolations defined in refs. |18|, M. 



In fact, from the point of view of gauge covariance the above definition satisfying (0) is not 
optimal. Instead of it one can also require the alternative condition 

f " dyA^y) = A Xfl . (9) 

J X 

This can be achieved, as one easily sees, by changing the definition of the continuation kernel 
in (|) to 

D,{x c - X ) = W ^e*^H^ , (10) 
il k 

where = 2sin(A; M /2). The consequence of eq. (§) is gauge covariance. Performing the gauge 
transformation on the lattice by the U(l) elements A x = exp(ia x ) and continuing a x to a 
function a(x c ) on the continuum in such a way that, for integer x c = x, 

a(x) = a x , (11) 

we obtain the relation for the gauge-transformed links 

UO? = KlfiU^K = exp S^ig dyA^\y)^ . (12) 

Here A^\x) is the gauge-transformed continuum field 

Al A \x) = A fl (x)-g-%a(x) . (13) 
A convenient interpolation satisfying eq. ([□]) is given by 

a{x c ) = ^D(x c - x)a x (14) 

X 

with the continuation kernel for scalar fields, analogous to in @ or (|10|), 

D(x c -x) = ^J2 eiHXc ~ X) ■ ( 15 ) 

k 



The gauge transformations become relevant if the effective action is not exactly gauge 
invariant (see the discussion in section |j). The choice of a gauge can be made on the basis of 
the minimality principle. This leads to the Landau gauge, which can be defined on the lattice 
by requiring 

0]^EE^ (is) 

X f_i=l 
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to be minimal with respect to gauge transformations. For U(l) gauge fields there are efficient 
algorithms to find this absolute minimum, for instance the one discussed in [[L{|, which will be 
used in the present paper. 

The advantage of the gauge-field interpolation given by eqs. (|5p, fllUp is its simplicity and the 
direct relation to momentum space, which will be useful for the evaluation of the determinant 
in momentum basis. Concerning topological charge, it defines a continuous gauge field on 
the torus which has a total classical topological charge zero. On large volumes this is not a 
serious constraint, because the parts of the volume can have any topological charge. In fact, 
the interpolated gauge field can be used to define a topological charge density operator with 
the appropriate renormalization procedure for composite operators. (See pOf and references 
therein.) 



3 Computing the determinant 

The Euclidean action for massless chiral fermions in the U(l) background gauge field A^x) is 
given in the continuum by 

S = Jd 2 x [$(x)j^iP(x) - igA^x^x)^ (PrQr + P l Ql) i>{x)} . (17) 

Here Pr = (1 + 73) /2 and = (1 — 73V2 are the chiral projectors for right-handed and 
left-handed fermions, respectively. We shall use the Pauli matrices for the 7-matrices in two 
dimensions: 7^ = cr^, (/1 = 1,2,3). Qr and Ql are the charges of the chiral fermion compo- 
nents. Note that in ( p|) we implicitly adopt the "doubling trick" [21, 22]: even if one of the 
charges Qr,l is zero, we use a Dirac fermion field. In this way the chiral fermion determinant 
is always a determinant indeed. 

The fermion matrix in momentum space corresponding to eq. (|17]) is 

Mll? L = ^<W7 • h - ig lfl (PrQr + PlQl) K-^ ■ (18) 
After multiplication by the fermion propagator we obtain 

A tQr,Ql — ajQr&l ' frO — a t^Qb,Ql / 1Q \ 

iV fc 2 jk! = M k 2 ki — — = d k2k! - & k2kl • iiyj 

This has the following matrix elements in spinor indices: 

iV Si Qi ( 1 ' X ) = <Wi ~~ QL(ak2-k u i - ia k2 ^ ku2 )(k ltl + ik lt2 )/kl , 

<f^(l,2) = 0, 
iVg^(2,l) = 0, 

Nglf L {2, 2) = 5 k2kl - Q R {a k ^ klA + ta k ,^ ku2 )(k hl - ik^)jk\ . (20) 

Here the explicit representation of the 7-matrices and the short notation a = gA/Q is used. 
For a massive vector-like fermion with mass m and Q = Qr = Ql, which will be used for 
Pauli- Villars fields, the matrix elements corresponding to ( p0|) are: 

iV^fci (1, 1) = h 2 ki ~ Q(a k2 -ki,i - iafca-fci^)(^i,i +ik lj2 )/{m 2 + k\) , 
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N k 2 ki (!> 2 ) = - im Q( a k 2 -k u i - ia k2 -k!,2) I '(m 2 + fc?) , 
( 2 > X ) = -imQ(a k2 - ku i + ia k2 ^ ku2 ) / (m 2 + k\) , 
Argf(2, 2) = 4 2fel - QK- fell i + ^ fc2 - fel , 2 )(A;i,i - ik 1)2 )/(m 2 + fc*) . (21) 

For the computation of the determinants of these matrices in momentum basis an appro- 
priate algorithm is the LU (lower-upper triangular) decomposition (see, for instance, [23]). It 
turned out to be both robust and sufficiently fast on the lattices considered. It can also be used 
for the computation of the full inverse matrix, and the algorithm can be organized in such a 
way that the matrix has to be stored only once. Of course, storing these large matrices even 
only once is the main limiting factor of the computation. For very large matrices also the time 
requirement is growing dangerously: it behaves as the third power of the matrix extension. 

In order to extend the range of feasible lattice sizes one can exploit some additional iterative 
procedures. Before describing them let us discuss the momentum cut-off scheme used. Since, 
according to the previous section, the Fourier components of the gauge field are constrained to 
the points (f|) of the Brillouin zone belonging to the gauge field lattice, it is natural to use a 
momentum cut-off for the calculation of the determinants of the infinite-dimensional matrices 
in (|19|)-([2T|). One can imagine to make the lattice finer for the fermions by adding more points 
to the gauge field lattice. In this case, however, the periodicity in momentum components is 
maintained, which introduces some non-zero elements also near the upper right and lower left 
corner, besides the ones near the main diagonal, for which a k2 - kx ,n ^ 0. This makes the effect 
of the cut-off stronger, therefore it is more advantageous to abandon periodicity and drop the 
extra non-zero elements. In this way, for momentum cut-offs much larger than ir (in gauge field 
lattice units), the matrix has a band structure. 

As a consequence of this band structure, one can effectively apply the iterative algorithm 
previously used for the numerical hopping parameter expansion in QCD ||24|| . For this one 



determines the traces of the powers of the hopping matrix K (here in momentum space). 
Having these traces one can use either the usual infinite expansion 



det(l-K)=exp -V^— , (22) 




or the finite polymer representation [E5| 

, , t ^ "^(-l)r^ 1 u ^+\ TrK^TrKP 2 Tr , . 

det(i -i^) = i + EEV- E ••• £ <W..+p,— — • • • — — • (23) 

u=lr=\ r[ pi=l p r =l Pi P2 Pr 

This latter is always convergent because it is a sum of a finite number of terms. In practical 
calculations one can go without any problems to j max ~ 100 in eq. (p2|) or to n max ~ 40 in 
eq. fl2"Bp. As a consequence of the band structure of the matrices K in (|iTj|)-(|2"TD, the storing of 



the full matrices is not necessary, and the computational load is growing as the second power 
of the matrix extensions times j max or n max . 

Inspection of the matrices in dl9|) -(^l|) shows that only the traces of even powers are non- 
zero. One can also easily see that 

Tr (K^^YX = Tr (k^A" . (24) 



5 



This corresponds to the relation 



det N Qr ' Ql ] * = det N Ql ' Qr . (25) 

Therefore the determinant of the vector-like fermion det is real. (One can also easily 

prove that it is non-negative.) 



Concerning the practical convergence of the trace expansions in (|22|) and/or (]23|) in the 
chiral Schwinger model with charges (3,4,5) the experience is negative. Typically neither of 
them converges, because in the calculated range the contributions rapidly increase. This is 
mainly the consequence of the large value of the charges (see next section). One can, however, 
easily save their advantages by calculating the determinant and inverse of the matrices iV 
truncated to a smaller sublattice (typically of the same size as the lattice for the gauge field), 
and then use 



det N = det iV small det [N~^ n (l - K)\ = det iV small det (1 - K new ) , (26) 

with 



K new = 1 - Aryi - K) . (27) 

It turns out that the expansions in the traces of the powers of Knew converge rapidly. The 
omitted higher-order terms could always be kept smaller than 10~ 8 , relative to the result. 



4 Convergence and gauge variance 

The commonly considered example of an anomaly-free chiral Schwinger model has fermion 
charges (Q R = 3,Q L = 0), (Q R = 4, Q L = 0) and (Q R = 0, Q L = 5). (The anomaly is 
cancelled if the sum of squared charges of the right-handed fermions is equal to those of the 
left-handed ones.) In order to regulate the ultraviolet divergence of the two-point function in 
a gauge-invariant way, one can introduce a Pauli-Villars vector-like fermion field with charge 
Q = 5 |26||. In this way the effective action S e g induced by the fermions is given by 



exp{-S eS (M PV )} = E cS (M PV ) = det N 3 ' det iV 4,0 det N 0,5 / det N 5{Mpv) . (28) 
Here Mpy is the mass of the Pauli-Villars regulator field in units of the gauge field lattice. As 



argued in ref. [14], in the continuum limit Mpy should be kept finite, for instance of order 1. 



This is necessary in order to maintain the possibility of a simple renormalization. 



Some insight into the behaviour of the effective action defined by eq. (p8|) can be obtained 
by numerically evaluating the determinants on some typical gauge configurations taken from 
Monte Carlo updating. For this I took quenched updating by the usual compact U(l) gauge 
field Wilson action. The gauge coupling has been fixed by (3 = g~ 2 = 8, which is a quite strong 
coupling for these fermions. Namely, the interaction strength is given by Qg, therefore weak 
couplings are beyond (3 ~ 25. The gauge configurations were transformed to Landau gauge 



by the algorithm described in [19]. The gauge field lattices were either 4 <g> 4 or 10 <S> 10. The 
lattices defining the momentum cut-off for the evaluation of the determinants always had an 
odd number of points (see discussion after eq. (§)), and they went up to 61 <S> 61. This means 
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that the momentum cut-offs went up roughly to 157r (in units of the gauge field lattice). For 
the Pauli-Villars mass values between Mpy = | and Mpy = n were tried. 

A first important question is how fast the infinite cut-off limit is reached by the effective 
action. It turned out that for every considered gauge configuration, taken randomly from the 
updating and transformed to Landau gauge, a good convergence could be achieved with the 
above cut-offs, provided that Mpy was not too large. For an illustration on 4 (g) 4 lattice see 
fig. [I]. The numerical results for nine configurations on 10 <8> 10 are shown in table [I]. 

It is interesting to investigate the gauge dependence of the determinants. The gauge trans- 
formation of the infinite-dimensional fermion matrix M in (|18D is given by A^MA, where in 
momentum space 

A fc2fcl = ^ / d'xe^- k2 > +ia ^ = 5 k2kl + ^&k 2 - kl + ^ E + • • • , (29) 



with &k denoting the Fourier components of ai(x). This infinite-dimensional unitary matrix is 
truncated by the momentum cut-off. Therefore, gauge invariance of the fermion determinant 
is lost even for vector-like fermions, which were gauge-invariant without truncation. The chiral 
fermion determinants remain gauge non-invariant also for infinite momentum cut-off. 
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fermion determinant on 4®4 

-i i i | i i i | i i i 



= 8 
M pv - 1.0 



(SO 



0.074 0.075 0.076 

Re(det) 



0.077 



Figure 1: The values of E e s(l) in the complex plane on a 4(g>4 gauge field with momentum 
cut-off on 11 g) 11 (triangle), 21 ® 21 (quadrangle), etc., up to 61 <g) 61 (eight-angle). 



The variation with gauge transformations is displayed in figure ^| for ten (3 = 8 gauge 
configurations on 10(8)10 lattice. The configurations were first transformed to Landau gauge and 
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ratios of fermion determinants on 10®10 



(8 = 8 



M pv -0.5 




0.5 1 1.5 2 

Rc( ratio) 

Figure 2: The residual gauge variance of the effective action shown by the (complex) 
ratios of _E e fj(|) atfer and before gauge transformation. The numbers are labeling different 
configurations. The circle with radius | around the point (1,0) is drawn to guide the eyes. 



then random gauge transformations were performed with parameters satisfying on the lattice 
points — 7r/20 < a x < tt/ 20. The gauge transformations were interpolated in the continuum 
as described in section |2|. The determinants were calculated on 31 ® 31 (momentum cut-off 
= 3.l7r). As is shown in the figure, the gauge variation is not very strong. In fact, both 
numerator and denominator of _E e fr(|) in fl2"%| ) always change by 4-5 orders of magnitude, but 
the ratio remains close to 1. Performing random gauge transformations with larger magnitude 
(up to — 7i < a x < 7r) shows an ever- increasing change in numerators and denominators, such 
that it becomes difficult to keep the numbers in the computer, but the overwhelming part of 
the variation is cancelled in the ratio. The cancellation can be further improved by taking more 
Pauli-Villars fields with appropriately chosen larger masses. 

There are two possibilities for dealing with this residual gauge variation of the effective 
action. First, one can try to tolerate it, keeping the momentum cut-off finite in gauge field 
lattice units. Second, more radically, one can enforce exact gauge invariance by defining the 
effective action to be equal to its value in Landau gauge along the whole gauge transformation 
orbit. The hope is that at the end, in the continuum limit, both these procedures lead to the 
same well defined theory. 

The gauge-field interpolation combined with momentum cut-off for the evaluation of the 
Pauli-Villars regulated determinants seems to work reasonably well in the (3,4,5) chiral Schwing- 



8 



er model. It can be expected that the effective action defined in this way leads to a well defined 
continuum limit. Of course, momentum cut-off is not the only possibility. Examples of other 
possibilities are, for instance, to take a finer lattice (in coordinate space) for fermions and to 



use the formalism of ref. [22] for the imaginary part of the effective action as suggested in [12 



or to take on the finer fermion lattice the SLAC derivative, as proposed by ref. [0]. One has 
to see which one of these (or some other) approaches has the most conceptual and practical 
advantages. 



Table 1: The values of the determinants on 10 ® 10 gauge field configurations 
with momentum cut-off 5. In. The first line for a given configuration label is 
the value on the "small" subspace det -/V small , the second line the correction 
factor obtained by trace expansion. For the configurations above the double 
line iVamau is with momentum cut-off 2.17T, below it with 3.17T. The complex 
numbers are given by pairs in parentheses. 





det N 3 ' 


det iV 4 >° 


det N 0,5 


det N 5 ^ 


1 


(2.1481,-0.8303) 
(1.0216,-0.0055) 


(0.3880,0.5783) 
(1.0364,-0.0288) 


(-68.1036,-40.6948) 
(1.0541,0.0308) 


3.5545 x 10 1U 
1.3258 


2 


(0.4173,-0.2723) 
(1.0374, -0.0053) 


(0.2104,-0.2211) 
(1.0672,-0.0148) 


(-0.0557,-0.1826) 
(1.1058,0.0367) 


5.5351 x 10 8 
1.4123 


3 


(0.0987,0.0139) 
(1.0482,-0.0000) 


(0.0169,0.0081) 
(1.0892,-0.0040) 


(0.0020,-0.0020) 
(1.1457,0.0155) 


5.6234 x 10 4 
1.4466 


4 


(0.6872,0.2887) 
(1.0303,0.0009) 


(1.2979,0.1028) 
(1.0561,0.0052) 


(1.2729,0.0435) 
(1.0929,-0.0166) 


9.3387 x 10 4 
1.3091 


5 


(0.3724,0.0593) 
(1.0268,-0.0073) 


(0.1504,-0.0241) 
(1.0458,-0.0117) 


(0.0073,0.0931) 
(1.0673,0.0132) 


1.3400 x 10 4 
1.2518 


6 


(0.2303,0.0165) 
(1.0271,-0.0008) 


(0.0764,0.0111) 
(1.0482, -0.0004) 


(0.0156,-0.0025) 
(1.0751,-0.0007) 


2.4780 x 10 1 
1.2412 


7 


(1.3350,-0.3054) 
(1.0474,0.0084) 


(2.0516,0.3677) 
(1.0884,0.0262) 


(-1.1129,-17.2413) 
(1.1481,-0.0658) 


3.6316 x 10 10 
1.5387 


8 


(0.00956,-0.00109) 
(1.0402,-0.0048) 


(0.000327, -0.000086) 
(1.0322,-0.6017) 


(0.0000037,0.0000019) 
(0.8859,-0.1220) 


5.920 x 10 8 
1.3736 


9 


(0.00945,0.00077) 
(1.0404, -0.0009) 


(0.000336,0.000054) 
(1.0257,-0.0070) 


(0.0000042, -0.0000024) 
(0.7122,0.2096) 


1.600 x 10 8 
1.3379 
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